Noise-Free Measurement of Harmonic Oscillators with Instantaneous Interactions 
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We present a method of measuring the quantum state of a harmonic oscillator through instanta- 
neous probe-system selective interactions of the Jaynes-Cummings type. We prove that this scheme 
is robust to general decoherence mechanisms, allowing the possibility of measuring fast-decaying 
systems in the weak-coupling regime. This method could be applied to different setups: motional 
states of trapped ions, microwave fields in cavity /circuit QED, and even intra-cavity optical fields. 
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Measuring the quantum state of a harmonic oscilla- 
tor, or, equivalently, its associated Wigner function 
is a fundamental task of quantum physics. Some propos- 
als allow a direct measurement, like propagating optical 
fields tested with homodyning techniques Q. Others, 
due to problems of accessibility, require indirect measure- 
ment schemes via interaction with a probe. This is the 
case of microwave fields in 3D cavities Q, circuit cavity 
QED with superconducting qubits [3], or the motion of 
trapped ions [5|. Different as they are, known techniques 
share a common problem: the noisy action of decoher- 
ence due to the probe-system finite interaction times. 

Recently, the reconstruction of a Wigner function in 
microwave cavity QED (CQED) was successfully real- 
ized [yjwith the aid of a dispersive probe-system inter- 
action [7| ■ Unfortunately, dispersive coupling is known to 
be slow and the required interaction time allows decoher- 
ence processes to disturb the measurement. In Ref. Q, 
a resonant method was proposed, but the Wigner recon- 
struction depends on the possibility of monitoring a few 
Rabi cycles, adding up to long probe observation periods. 
In the case of trapped ions, the measurement techniques 
are quite similar and a recent experiment Q made use of 
numerical integration over several Rabi cycles to achieve 
the goal. These long interaction times are particularly 
harmful in the case of fast-decaying systems. For exam- 
ple, state reconstruction of intracavity optical fields has 
not been experimentally attempted, to our knowledge, 
due to their weak coupling with atomic probes. 

In this paper, we present a method to measure the 
quantum state of a harmonic oscillator through instanta- 
neous probe-system interactions 13, preventing de- 
coherence from disturbing the measurement. The har- 
monic oscillator is allowed to interact with a two-level 
probe for an arbitrarily short time via a selective interac- 
tion 12], 13, EH in the Jaynes-Cummings (JC) model [l5j |. 
The information is then collected from the second time 
derivative of the probe population at zero interaction 
time. The scheme permits to measure the population 
field distribution and, with the support of coherent dis- 



placements, the associated Wigner and Q functions at 
any point in phase space with arbitrarily small influence 
of decoherence. From this data the full Wigner func- 
tion can be reconstructed by a simple fit as in [6j or by 
more sophisticated techniques [3], e.g., involving maxi- 
mum likelihood estimation (l7| and taking into account 
the imperfections of the measurement process [2j]. 

Typically, a selective interaction can be built when a 
three-level probe, driven by a classical and a quantized 
field, is reduced to two metastable states after adiabatic 
elimination of the third level, allowed by a large detun- 
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ing A. The associated Hamiltonian reads 



H eS = h^\g)(g\ 
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+H 



(\9)(e 



\e)(g\a) 



(1) 



Here, {\g), |e)} are the (metastable) ground and excited 
states of the two- level probe, {a, a'} are the harmonic os- 
cillator annihilation and creation operators, respectively, 
and = Slif^/A is the effective JC coupling strength. 
The first and second terms on the r.h.s. of Eq. ([1]) are 
AC-Stark shifts associated with each of the probe levels, 
the second one depending on the number of oscillator 
excitations. This means that any effort at tuning the 
JC coupling to resonance will succeed only for a selected 
JC doublet H N ■ {\g)\N), \e)\N - 1)}, leaving all other 
doublets, for which n ^ N, slightly or completely off- 
resonance. It can be shown that, under proper tuning of 
the excitation fields, the condition CI2 3> fli^/N assures 
neat selectivity in the JC model, where resonant Rabi 
oscillations will happen only inside the subspace TLn . In 
this case, Eq. ([1]) turns into the selective Hamiltonian 

H N = hVNn(\g)(e\ ® \N)(N - 1| + \e)(g\ ® \N - 1)(N\), 



(2) 



describing the flip-flop interaction of two effective spin- 
1/2 systems, {\g), |e)} and {\N-1), \N)}. Specific imple- 
mentations have been proposed in microwave CQED |l2| 
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and in trapped ions [l3(, but other systems, like an atom 
inside an optical cavity, or a superconducting qubit cou- 
pled to a coplanar waveguide resonator, can also enjoy a 
similar behavior. 

From the unitary time evolution of the total density op- 
erator, p = [HN,p]/ih, the first and second time deriva- 
tives of the expectation value of a time-independent 
probe operator B can be expressed as 



d{B) 1 i\r h n 



d 2 (B) 
dt 2 



1 



■{[[B,H N ],H N ]), 



(3) 
(4) 



We study here a more general case, allowing decoherence 
in the (field) system but disregarding decoherence in the 
probe. We do that based on the fact that most physical 
setups use probes with long lifetimes compared to the 
ones of the systems to measure. Under this assumption, 
we consider the mo st g eneral master equation p = Cp in 
the Lindblad form [181 ], 
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+ /J Wm (AnMm — ^Al n A m p — -pA m A m ) , (5) 



where n m express decay rates and Lindblad operators A r 
and Ajn are associated with the (field) system. 
We calculate first 

d{B) — ^ (\B H ]\ 
~~dt~ ~ih 



+Tr 



Km (A m pA m — — A m A m p ——pA m A m )B 



The second term on the r.h.s. can be rewritten under the 
trace as 



1 



Tr- 



,B]A m + A m [B,A r , 



and, since [-B,^J = [B, A m ] = 0, we have 



dt in 



(6) 



In this expression, see similarity with Eq. ([3]), the dy- 
namics of probe operator expectation value (B) does not 
seem affected by the system decay. This is certainly not 
the the time-dependent expectation value on the 

r.h.s. of Eq. ([6]) will in general be susceptible to deco- 
herence: the calculation involves time-dependent p(t) fol- 
lowing Eq. ([5|). However, at t — the noise terms vanish 
identically, i.e., the time derivative of any probe expecta- 
tion value at t = is independent of any field-decohering 



Lindblad environment. For this particular time, choos- 
ing a probe operator B = |e)(e| and probe-system initial 
state p(0) = \g)(g\ ® pj, we obtain 



dP e (r) 



dr 



0, 



T = 



with dimcnsionless time t = VLt and Pc{t) = ( |e) (e| 
the second derivative, we obtain 



(7) 
For 



^^=j^{[[B,H N ],H N ])+Tr Y,K m (A m f>Al 

-\A m A m p-^pA m A m )[B,H N ] . (8) 

This is a different situation and the trace will not van- 
ish in general, as before. However, considering again 
the observable B — |e)(e| and probe-system initial state 
p(0) = \g)(g\ <8 Pf, it follows that the second term in 
Eq. ||5J) vanishes again, and we have 



d 2 P e (T) 



dr 2 



= P, 



N • 



(9) 



T=0 



Here, Pm — T^\pf\N)(N\\ is the probability of finding 
Fock state \N) in the initially unknown harmonic oscil- 
lator state. 

Equation ([9]) describes a remarkable result, it shows 
that the curvature of the function P e (r), at vanishing 
t = 0, contains undisturbed information about P/v- This 
valuable field information is encoded correctly in the level 
statistics of the two-level probe even in presence of a field 
reservoir of a general kind. Needless to say, all previous 
results hold when the usual thermal bath is considered 
as a reservoir for the harmonic oscillator. At zero tem- 
perature, for example, Lindblad operators A m and JV m 
would have to be replaced by a and , respectively, and 
k would represent the decay rate of the single field mode. 
The counter-intuitive results of ([6]) and © are of an in- 
finitesimal nature and, without harming their theoreti- 
cal importance, should suffer high-order corrections when 
dealing with a discrete sampling of interaction times, as 
will be explained later. 

In order to measure the complete field population, P n , 
Vn, one just needs to tune resonantly the other selected 
subspaces Tl n and follow a similar procedure. The mea- 
surement of all P n allows the estimation of the complete 
Wigner function W(a) of p, conditioned to the realiza- 
tion of previous arbitrary field displacements D(—a) in 
phase space. For that, we have to recall that the Wigner 
function can be expressed among other possibilities, 
as 



W(a) 



n=0 



l) n P n (-a), 



(10) 



where P n (—a) stands for the field population after a dis- 
placement D(—a). Another rather original manner of 
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reconstructing the quantum field state in phase space is 
via the instantaneous measurement of the Q- function [l[ , 
defined as Q = (a\pf\a). It can be shown that 

Q(a) =Tr[D(-a)p f D\-a)\0)(0\] = Tr[ P/ (- a )|0)(0|] 

= P (-a). (11) 

This means that, following Eq. measuring instanta- 
neously the probability of having Fock state |0) after field 
displacements D(—a), amounts to a full measurement of 
the Q-function. Note that this will only require the tun- 
ing of a single selective subspace, reducing enormously 
the experimental efforts when compared to the case of 
the Wigner function. Typically, coherent displacements 
can be realized at a very high rate, depending mainly 
on the intensity of the excitation fields, so this is not a 
critical issue. 

From a fundamental point of view, our proposal sug- 
gests that, no matter how short the lifetime of a certain 
system is, there would always exist the possibility of en- 
coding its quantum information in two-level probe statis- 
tics at infinitesimal interaction times. In other words, 
measuring quantum states may not require long-living 
systems or strong probe-system coupling, against com- 
mon belief. In this work, we have not proved this con- 
jecture in general. However, we have given a particular 
example, the case of a quantum harmonic oscillator or 
a single mode field, that is applicable to many physical 
setups. 

Turning to more practical considerations, we stress 
that our scheme employs exclusively, as a final readout 
mechanism, the measurement of the population of the 
excited state of a two-level probe at different times. This 
probe population is measured directly by ionization, in 
the case of cavity QED, or fluorescence, for trapped ions 
and atoms, with suitable additional fields, different from 
the ones used in Eq. (1). This is a standard measurement 
in many quantum-optical experimental setups of interest 
and has been realized routinely with high accuracy and 
efficiency (EHGit 

Estimating derivatives from a finite-time experimental 
sampling will produce higher order corrections and an 
overhead in the number of repetition measurements. In 
consequence, the benefit of instantaneous measurements 
demands an improved measurement accuracy, which can 
be seen as follows. To determine the discrete second 
derivative in Eq. ©, P e (0), we measure P e at 0, r, 2r, 
for small t, and calculate 



P e (2r)-2P e (r)+P e (0) 



P e (0) + o(r), (12) 



The measurement results are scattered around each P e 



with variance A 2 



A 2 , referring to quantum- 



mechanical and technical noise, the latter arising from 
imperfections in preparation, timing, and measurement. 
For small values of r and probe operator B — |e)(e|, we 



have A 2 (r) = (B 2 ) - (B) 2 = P e (r)-P e (r) 2 ~ P N r\ fol- 
lowing the Taylor expansion of P e (r) around zero time. 
Performing M measurements at each of the three times 
of Eq. (|12|) . we can reduce the uncertainty in the three ex- 
pectation values by a factor l/yM. Thus, our estimate 
for the second derivative comes with an error variance 

. 6Pvr 2 + 4A 2 



A 



Mr 4 



(13) 



while the signal is ~ Ppj. Then, in order to achieve a 
signal-to-noise ratio larger than unity, P/v/A > 1, wc 
need a number of measurements 



M > T~ 



(14) 



approximately. In this way, as we reduce r by a factor /, 
in order to improve the approximation of P e (0) , we need 
to increase the number of measurements by / 2 (or / 4 if 
technical errors dominate) to maintain the desired signal- 
to-noise ratio. Consequently, the total probe-system in- 
teraction time Mt summed over all measurements in- 
creases, while each measurement outcome represents the 
effect of an arbitrarily short interaction time r. Although 
this is a high price to pay, contamination of the oscillator 
by decoherence can be kept arbitrarily small, no matter 
how fast it is. 

Quantitatively, the determination of P e (0) is associ- 
ated with an error ~ P e (0)r due to the finiteness of r. 
In this case, we can use Eq. (J5j to obtain an expres- 
sion for P e (0) and observe two contributions: a unitary 
component, present even in the absence of noise, and an 
additional term cx k/Q due to the noise [2(ij ]. In order 
to determine P e (0) to accuracy q, we need r < q/P e (0). 
This, in turn, requires M > q~ 2 P e (0) 2 for ideal measure- 
ments and M > g~ 4 P e (0) 4 A 2 in the case of technical- 
dominated errors. Hence, the number M of required 
measurements increases polynomially with k (in the case 
k 3> Q like k 2 resp. k 4 ). Remarkably, no matter how 
strong k is, we can always make sure that it does not 
affect the measurement result by shortening r. This pro- 
cedure finds its natural limit r > Q/uj, where cj is the 
smallest dominant frequency of the specified dynamics 
(u> = fi|/A in our example, which implies r > f^/f^)- 
Note that we only need to know an upper bound of the 
decoherence rate to perform an accurate measurement, 
while alternative methods involve specific decoherence 
models. 

The proposed method can be applied to any physical 
system enjoying JC selective interactions, and may allow 
for the complete measurement of elusive fast-decaying 
systems, as is the case of intracavity optical fields, among 
others. Furthermore, since the scheme is based on in- 
finitesimal transfer of information, this method could ob- 
serve initial state conditions. For example, it could be 
used to test the initial purity of a system, as well as its 
unavoidable entanglement with external degrees of free- 
dom at zero interaction time. This may prove useful for 
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testing the validity of usual assumptions taken in deco- 
herence theories, such as system and environment being 
initially in a separable state or markovian dynamics. 
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